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ABSTRACT 


The spacecraft in the next series of Geostationary Operational Environmental 
Satellites (GOES-Next) are Earth pointing and have 5-year mission lifetimes. Be- 
cause gyros can be depended on only for a few years of continuous use, they will 
be turned off during routine operations. This means attitude must, at times, be 
determined without benefit of gyros and, often, using only Earth sensor data. To 
minimize the interruption caused by dumping angular momentum, these spacecraft 
have been designed to reduce the environmental torque acting on them and incor- 
porate an adjustable solar trim tab for fine adjustment. A new support require- 
ment for GOES-Next is that of setting the solar trim tab. Optimizing its setting 
requires an estimate of the unbalanced torque on the spacecraft. These two re- 
quirements, determining attitude without gyros and estimating the external torque, 
are addressed by replacing or supplementing the gyro propagation with a dynamic 
one, that is, one that integrates the rigid body equations of motion. By process- 
ing quarter-orbit or longer batches, this approach takes advantage of roll -yaw 
coupling to observe attitude completely without Sun sensor data. Telemetered 
momentum wheel speeds are used as observations of the unbalanced external torques. 
GOES-Next provides a unique opportunity to study dynamic attitude propagation. 

The geosynchronous altitude and adjustable trim tab minimize the external torque 
and its uncertainty, making long-term dynamic propagation feasible. This paper 
presents the equations for dynamic propagation, an analysis of the environmental 
torques, and an estimate of the accuracies obtainable with the proposed method. 
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1 . INTRODUCTION 


Accurate attitude determination typically requires a large amount of data taken at 
different times. This takes advantage of averaging to reduce the effects of sensor 
noise but requires a means of attitude propagation. Three-axis stabilized space- 
craft usually carry gyros that measure how much the spacecraft rotates over short 
time intervals. For the next series of Geostationary Operational Environmental 
Satellites (GOES-Next), however, the gyros will be turned off when the spacecraft 
is on-station, making the usual method of attitude determination impossible. 

GOES-Next has Earth and Sun sensors, and when both provide data, it is possible to 
compute a "single-frame" attitude solution. The Earth sensors provide pitch and 
roll data; the Sun sensors provide pitch and yaw data. However, the Sun is visi- 
ble to the Sun sensors for only two-thirds of the day-long orbit, causing an 8-hour 
period each day when yaw cannot be observed directly. Having an alternative to 
gyro propagation would make it possible to compute the yaw when the Sun is visible 
and then predict it for later times when the Sun is out of view. The obvious can- 
didate for this role is the dynamic equation for rigid body rotation, or Euler's 
equation. 

Using Euler's equation to propagate for attitude estimation is not a new idea. In 
1976, Lefferts and Markley (Reference 1) and Markley and Wood (Reference 2) applied 
dynamic propagation to Nimbus-6. The estimator included detailed torque and dy- 
namic modeling and, in addition to solving for attitude and angular velocity, al- 
lowed the estimation of torque model parameters. It worked well with simulated 
data but was unable to duplicate real Nimbus-6 attitude histories. This difficulty 
seemed due to imperfectly known environmental torques, unmodeled control system 
activity, and uncertainty about the mass properties of the spacecraft. 

Fein (Reference 3) concentrated on the Idea of estimating environmental torques 
from wheel speeds and sensor data. He attributed short-term variations in the 
speed to the control system and long-term variations to the environmental torques. 
Based on the knowledge that the attitude remained close to nominal, he was able to 
model the torques with low-order polynomial functions. Sensor observations served 
to correct the propagated attitude. Although only time spans up to 22 minutes 
were considered, agreement with observed attitude histories was good. 
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Because GOES-Next is at the much higher geosynchronous altitude, environmental 
torques are expected to be much smaller and correspondingly less uncertain. This 
should make dynamic propagation more feasible for GOES-Next than for Nimbus-6. 
Prospects for GOES-Next are further encouraged by its large pitch momentum bias, 
which stabilizes the spacecraft yaw and roll by keeping the pitch axis aligned 
with the orbit normal (Reference 4). 

This paper describes an attitude estimator that uses Euler's equation, adapts it 
to the GOES-Next mission, and estimates its accuracy for yaw determination. 

2. GOES-NEXT ATTITUDE 


The attitude of GOES-Next is defined relative to a rotating reference coordinate 
system with its z-axis pointing to the center of the Earth; y-axis in the direction 
of the negative orbit normal; and x-axis oriented so that x, y, and z form a right- 
handed orthogonal triad. The spacecraft attitude relative to the reference system 
is defined by a 3-1-2 Euler axis sequence with the three rotation angles referred 
to as yaw (y), roll (r), and pitch (p) (Reference 5). 

The spacecraft roll, pitch, and yaw axes (x, y, and z, respectively) are close to 
the principal axes of the spacecraft moment-of-inertia tensor. The diagonal compo- 
nents of this tensor in the roll, pitch, yaw frame are I y = 3364.376 kilogram- 
meters squared (kg-m 2 ), I = 954.936 kg-m 2 , and I z = 3461.393 kg-m 2 , and all of 
the off-diagonal elements (the products of inertia) are less than 30 kg-m 2 in 
magni tude. 

Estimated attitude uncertainty is due to sensor noise, sensor visibility, and prop- 
agation noise. The standard deviation of the sensor noise is a measure of the 
accuracy of the sensors. Visibility, in this case, involves only Sun visibility 
since the Earth is assumed to be in view at all times. Propagation noise includes 
errors in the models of environmental and control torques and unmodeled contribu- 
tions to these torques. 

2.1 SENSORS AND SUN VISIBILITY 

GOES-Next has Earth sensors that measure pitch and roll and Sun sensors that meas- 
ure pitch and yaw; these measurements are telemetered at 0.512-second intervals 
(Reference 6). Under normal conditions, the Earth sensors provide continuous 
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measurements. The Sun sensors, however, do not cover the spacecraft z-axis and so 
do not see the Sun around local midnight. 

There are two Earth sensors that scan east-west across the Earth disk and combine 
their measurements to give pitch and roll.. Because of the high altitude, horizon 
height uncertainty contributes little to sensor error, and the standard deviation 
of these observations is given as 0.01 degree (la). The quantization is also 
0.01 degree. 

There are two types of Sun sensor. The digital Sun sensor (DSS) is the more accu- 
rate of the two, having standard deviation of 0.0042 degree Oct) and quantiza- 
tion of 0.125 degree. There are also two coarse analog Sun sensor (CASS) systems. 
These are less accurate, having unspecified noise characteristics and quantization 
of 0.039 degree. The error for these sensors is on the order of 1 degree on the 
boresight and increases away from the boresight. 


GOES-Next rotates once a day about its y-axis, which is roughly 67 degrees from 
the ecliptic plane. The percentage of each orbit having Sun coverage is approxi- 
mately the fraction of the spacecraft x-z plane in the Sun sensor fields of view. 
As seen from Figure 1, this leaves about one-third of each day without coverage 


SUN NOT 
VISIBLE 



Figure 1. Sun Sensor Field of View 
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and without yaw observations. It is the more accurate DSS that sees the Sun just 
before the loss of coverage. 

2.2 TORQUES 

The torques acting on a spacecraft arise both naturally and from the control sys- 
tem. Natural disturbances include solar radiation, gravity gradient, residual 
magnetic dipole, and aerodynamic torques. For GOES-Next, which is at geosynchro- 
nous altitude, atmospheric torques are negligible. The torques caused by the con- 
trol system include those due to the magnetic torquers, changing wheel speeds, and 
thruster activity. Because the thrusters are not expected to be used more than 
once a day, they are not modeled here. 

Solar radiation causes the largest environmental torques on GOES-Next (Refer- 
ence 7). There is a controllable flap at the end of the solar array that is com- 
manded from the ground to minimize the solar torque. Proper setting of the trim 
tab can reduce the solar torque to 10“ 7 newton meters (N«m), but the residual 

_5 

torques may be as large as 10 N*m. The solar torque is approximately a constant 
scalar times the cross product of the vector from the center of pressure of the 
solar array to the center of pressure of the solar sail, which is closely aligned 
with the spacecraft pitch axis, and the Sun-to-spacecraft vector. Thus, the solar 
torque vector is nearly constant in inertial space and is mostly in the spacecraft 
roll/yaw plane. In the spacecraft body frame, the largest components of the solar 
torque are the roll and yaw components, which have sinusoidal time dependence with 
the orbit period (one sidereal day) with nearly equal magnitudes, and a 90-degree 
difference in phase. Figure 2 illustrates this dependence of the roll and yaw 
torques. 

Gravity gradient torques depend only on the spacecraft attitude and moment-of- 
inertia tensor, so they are easily modeled. The largest errors in the gravity 
gradient torque model come from uncertainty about the inertia tensor for the 
spacecraft. 

Magnetic torques come from residual and control magnetic dipoles (References 8 and 
9). The residual dipoles are due to electric currents, magnetized material in the 
spacecraft, and commanded torquer biases. They are considered to be almost con- 
stant in the spacecraft body frame and to have magnitudes of about 8 ampere-meters 
2 

squared (A-m ). At geosynchronous altitude, this can produce a torque of at most 
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Figure 2. Solar Radiation Torques 

10 -7 N«m. Because of the equatorial orbit of GOES-Next, the Earth's magnetic field 
is mainly along its pitch axis. Therefore, the torque due to residual dipoles is 
approximately constant in the body frame. The control dipoles vary slowly over 
time, producing torques as large as 10 -6 N*m. 

2.3 MOMENTUM WHEELS 

GOES-Next is equipped with two momentum wheels with their axes in the pitch-yaw 
plane, canted at an angle of 1.656 degrees from the pitch axes. These produce a 
pitch angular momentum bias of 

H = -I («. + « 0 ) cos (1.656°) (1) 

W 1 2 

and a yaw angular momentum of 

h = I («. - «,) sin (1.656°) (2) 

W 1 2 
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where w, and w 9 are the angular velocities of the two wheels in radians per second 

I c 2 

(rad/sec) and I w - 0.1082 kg-m is the moment of inertia of each momentum wheel. 

The nominal on-orbit wheel speeds are 5485 revolutions per minute (rpm) for both 
wheels, giving H - -124.2 Nms and h - 0. Commanding the wheels in the same direc- 
tion gives a pitch control torque, and commanding them in opposite directions gives 

a yaw control torque. The dominant error contribution of the momentum wheels is 

-7 2 

the torque ripple, which has a maximum spectral density of 6 x 10 (N«m) /hertz 

(Hz) for each wheel (Reference 10). The momentum wheel speed is measured by a 
tachometer, sampled at 0.512-second intervals, with a quantization of 0.0163 rpm 
and a noise level of 0.0978 rpm (3<x) (Reference 11). 

Should one momentum wheel fail, the other wheel is operated with a smaller reac- 

2 

tion wheel (moment of inertia = 0.008626 kg-m ) whose axis is along the yaw axis. 
The speeds of these wheels are nominally set to give a smaller pitch angular mo- 
mentum bias and zero net angular momentum along the yaw axis, as before. The 
present analysis is easily extended to this backup situation, but it will not be 
considered further. 


3. DYNAMIC MODEL 

For an Earth-pointing spacecraft in a circular orbit, Euler's equation and the 
kinematic attitude equations can be expressed in terms of roll, pitch, and yaw; 
the resulting equations can be linearized for small deviations from nominal atti- 
tude (e.g., pages 608-610 of Reference 5). The resulting equations for GOES-Next, 
taking into account the rotation of the reference frame (at one revolution per 

_c 

sidereal day, or co Q = 7.29 x 10 rad/sec) are 

I x r - (o Q [H - 4(1 - I z ) « Q ] r - CH + (I x - I y + I z ) « Q ] y - N x + o> 0 h (3a) 

IyP + 3»^(I X - I z > p - N y - H (3b) 

I z y - U() [H - <I y - I x ) « 0 ] y + CH + (I x - I y + I z > « Q ] ? « N z - N w (3c) 
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where r, p, and y are the spacecraft roll, pitch, and yaw in radians; a dot de- 
notes a time derivative; I , I , and I are the diagonal components of the space- 

x y 2 

craft moment-of-inertia tensor (the off-diagonal elements have been neglected); H 

and h are the internal pitch and yaw angular momentum (as defined in Section 2.3); 

N , N , and N are the external torques on the spacecraft (excluding the gravity 
x y z 

gradient torque, which is included on the left-hand sides of the equations); and 
N w is the yaw component of the momentum wheel control torque, so that 

h = N w (4) 

Equations (3a), (3b), and (3c) exhibit the well-known fact that the spacecraft 
pitch motion is independent of the coupled roll/yaw motion, to within the accuracy 
of the linear approximations used to derive these equations. The GOES-Next pitch 
is well determined by continuous Earth sensor measurements, as discussed in Sec- 
tion 2.1, so this paper will concentrate on the roll /yaw motion. 

The numerical values for I and I z given in Section 2 establish the validity 
of the approximation 

I x - I z - I = d x I z ) 1/2 = 3412.54 kg-m 2 (5) 

The differences of the principal moments of inertia can be ignored in the coeffi- 
cients containing the bias momentum in Equations (3a) and (3c) since 

“ 0 dy - I) - 0.0015 H (6) 

With these approximations. Equations (3a) and (3c) can be written 

r * u 0 « n r + (o> n + w 0 ) y + (N x + u 0 h)/I (7a) 

y - co 0 co n y - U n + ( 0 o > r + (N z - N w )/I (7b) 

where the nutation frequency « n is given by 

o) = H/I = 0.0364 rad/sec (8) 

n 
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This corresponds to a nutation period of 173 seconds. 

Precise attitude estimation will require simultaneous estimation of unmodeled 
torques, in view of the difficulty of accurately modeling the disturbance and con- 
trol torques. As discussed in Section 2.2, the principal torque modeling errors 
have both a constant term N c and a sinusoidal term N p with angular frequency 
(ji Q . Thus, the torques are 


1 N„ - 

N* + 

+ 

N^„ + w„ 

(9a) 

1 X 

X cx 


px X 


N = 

N* + N 

+ 

N _ + w_ 

(9b) 

z 

z cz 


N 

N 

Q. 


1 

N = N* 
w w 

+ 

w w 

(9c) 


where n£, n£, and n£ are the modeled roll, yaw, and wheel torques, respec- 
tively, and w , w , and w are independent white noise processes. The estimated 

A ^ W 

values of N c and N p will include all the unmodeled torques with the prescribed 
time dependence. 

A nine-component state vector containing all the parameters to be estimated is 


x = [r, r/« n , y, y/« n , h, N px , N pz , N cx , N CZ ] T (10) 

where superscript T denotes the matrix transpose. For GOES-Next, variations in H, 
and therefore in « n , are very small and can be ignored. Then, combining Equa- 
tions (4), (7), (9), and (10) with the assumed sinusoidal dependence of N p gives 
the dynamic equation 


x«Fx+u+w 


where 


F 


A J B 

— | — 
o s c 


(11a) 


(lib) 
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with 


0 0 

“o 0 0 

0 0 0 

0 -(c o + w ) w 

n 0 0 

0 0 0 0 

-O 1 0 1 

0 0 0 0 

0 0 10 




and 0 is a 5-by-4 matrix of zeros. 

The nine-component vectors u and w are given by 


u - [0, NJ/H, 0, (NJ - N*)/H, NJ. 


0 

Gd + W_ 

n o 

(l) 

n 

0 


0 

0 

0 

1 


0 

0 

0 

0 

0 


0, 0, 0, 0] T 


(11c) 


(lid) 


(lie) 


(12a) 


and 


w - CO, w x /H, 0, (w z - w w )/H, w w , w px , w pz , w cx , w cz ] T (12b) 

where w px , w pz , w cx , and w cz are independent (of each other and of w x , w z , and w w ) 
random noise processes. If w and w are not zero, N will not be truly constant 

C X CZ v 

but will behave as a random walk. Nonzero w and w will give random walk de- 
pendence to the phase and amplitude of N p . 
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Equation (11a) has the formal solution 


x(t) 


4>(t, t n ) x( t ) + 
0 0 



<t>(t, t')[u(t') + w(t')] dt* 


where the state transition matrix 4*< t , t 0 ) is a solution of 

*<t, t Q ) = F *(t. t Q ) 


(13) 


(14a) 


with the initial condition 

4>(t Q , t Q ) = I g = the 9x9 identity matrix (14b) 

To avoid confusion, capital I denoting an identity matrix will always have a 
numerical subscript, and capital I denoting moment of inertia will have either a 
literal subscript or no subscript. The state transition matrix has a partitioning 
similar to Equation (lib): 


*(t. v 


♦(t. t 0 ) ; e(t, t 0 ) 
o v(t, t 0 ) 


(15) 


Substituting Equation (15) into Equations (14a) and (14b) gives 


4»(t, t Q ) - A *(t, t Q ) with 4>(t 0 , t Q ) - I 4 


(16a) 


t(t. t n ) - C t(t. t n ) with ^(t 0 , t Q ) = I 


(16b) 


and 


G(t, t 0 ) = A 6(t, t Q ) + B *(t, t Q ) with 0(t o , t 0 ) - 0 


(16c) 
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The solution of Equation (16c) is 


ect, t 0 ) 


J *(t, t- 


) B f(t\ t Q ) df 


so 0(t, t Q ) is determined when <J>( t , t Q ) and f(t, t Q ) are determined, 
lution of Equation (16b) is easily seen to be 


Y<t, t 0 ) - 


1 

0 

0 

0 

0 


0 

C c 

- S C 

0 

0 


0 

s„ 


0 

0 


0 

0 

0 

1 

0 


0 

0 

0 

0 

1 


where 


c = cos « (t - t ) 
o o o 


and 


s = sin to (t - t ) 
o oo 

It is thus seen that the structure of the matrix C gives the desired 
dependence of N p . 

Equation (16a) is more difficult to solve. Since A is constant. 


4>(t, t Q ) - exp [A(t - t Q )] 


The characteristic equation of A is 


0 - det [A - XI 4 ] - | X 2 + to p j |x 2 + to 2 j 


(17) 

The so- 

(18) 

(19a) 

(19b) 

sinusoidal 

( 20 ) 

( 21 ) 
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so the eigenvalues of this matrix are ± i« Q and ± i« n . With Equation (20), this 
shows that <|>(t, t Q ) comprises sinusoidal terms with angular frequencies co Q and 
w n - This is the justification for referring to « n as the nutation frequency. 

Note that the characteristic equation is intractable without the assumptions lead- 
ing from Equations (3a) and (3c) to Equations (7a) and (7b). With the periodic 
nature of <|>(t, t Q ) known, it is not too difficult to show that the solution of 
Equation (16a) is 


♦<t. y - yt, y ♦ yt. y 


(22a) 


where 


*o (t ’ t o ) 


CO - <0 

n o 


♦n (t - V 


co - co_ 
n o 


with 


and 


CO C 

n o 

-co S _ 

n o 

o 

00 

c 

3 

coC 

n o 

— coS 
0 0 

— coC 
0 0 

co C 
0 0 

-uS 
0 0 

— co S 
n o 

—to C 

no 

<oC 

n o 

-co 

n o 

-coC. 
_ 0 0 

coS 
0 0 

-co S 
0 0 

-co. c 

o o_ 

ions (19a) and 

(19b), and 


— coC 

o n 

<*) s 
n n 

-cos 

o n 

— co C 
n n 

OdS 

o n 

n n 

-<o.C_ 
o n 

co S 
n n 

to S_ 

o n 

w_c 
n n 

-co c„ 

o n 

co_S_ 

n n 

_"o c „ 

—co s 
n n 

co S„ 

o n 

coC 
n n 

c n = cos 

<0 ( t 
n 

- V 


s n 5 sin 

CO (t 

n v 

- V 



(22b) 


(22c) 


(23a) 


(23b) 
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Substituting Equations (18) and (22a), (22b), and (22c) into Equation (17) gives 
integrals that are convolutions of trigonometric functions and can be evaluated 
without too much difficulty. The terms in 0(t, t Q ) have only six kinds of time 
dependence: l-c n , s n , l-c Q , s Q , « 0 (t - t Q )c 0 , and » Q (t - t Q )s 0 . 

4. OBSERVABILITY 

Before proceeding to detailed treatment of the attitude estimation algorithms, it 
is useful to establish that the attitude state defined in Section 3 is observable. 
To discuss observability in a more general sense, consider an n-dimensional state 
vector and an m-dimensional vector of measurements of x: 


g - G x (24) 

where G is an m-by-n matrix. Then the state is observable if and only if the 
nm-by-n observability matrix 


M = 


G 

GF 

GF 2 


GF 


n-1 


(25) 


has full rank n, where F is the dynamic matrix defined by Equation (11a) for the 
n-component state x. 

Consider first the full nine-component state, with F given by Equations (lib) 
through (lie). If roll, yaw, and wheel tachometer measurements are available, 
m - 3 and 



0 0 0 
0 1 0 
0 0 0 


0 0 0 0 0 
0 0 0 0 0 
1 0 0 0 0 


( 26 ) 
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A computation of H shows that the state is observable in this case; in fact, com- 
puting G, GF, GF , and GF achieves full rank. Thus, the roll/yaw attitude, atti- 
tude rates, wheel speed, and torques are all observable with roll, yaw, and 
tachometer measurements. 

The observability when yaw measurements are not available is also of interest, as 
discussed in Section 2.1. With only roll and tachometer measurements, m = 2 and 


G" = 


1 

.0 


0 0 0 
0 0 0 


0 0 
1 0 



(27) 


The observability matrix computed with this has rank eight, so the nine-component 
state is not observable. 

To ascertain which parameters are observable and which are not, consider the eight- 
component state vector obtained by deleting the constant yaw torque from xf: 


x 1 = [r, r/u n , y, y/to n , h, Np X , Np Z , N CX ]T 
This obeys the state equation 


where 


with 


x‘ = F' x 1 + u‘ + w 1 



0 


B' 


1 

H 




0 

0 


0 0 
1 0 

0 0 
0 1 


0 

1 

0 

0 


(28) 


(29a) 


(29b) 


(29c) 
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and 


C* 



0 0 

0 1 

-1 0 

0 0 


0 

0 

0 

0 


With only roll and tachometer measurements 


G' 


1 0 0 0 0 0 0 
0 0 0 0 10 0 


0 

0 


(29(1) 


(30) 


Forming the observability matrix from products G'(F') for power k up to six gives 
full rank, establishing that the attitude, attitude rates, wheel speed, periodic 
torque and constant roll tofque are observable with only roll and tachometer meas- 
urements. It is more difficult to observe the reduced state without yaw measure- 
ments than to observe the full state using yaw measurements along with the roll 
and tachometer measurements, as indicated by the need for higher powers of F in 
the former case. 

The preceding computations establish that the full nine-component state is observa- 
ble with roll, yaw, and tachometer measurements, but the constant yaw torque is 
unobservable in the absence of yaw measurements, the other eight components of the 
state remaining observable from roll and tachometer measurements. This indicates 
that the constant yaw torque must be estimated during periods with yaw measure- 
ments, and that the errors in this torque must not grow too rapidly in periods 
without Sun sensor visibility. This is equivalent to the assumption that the error 
source w cz is small . 


5. KALMAN FILTER 

In principle, the GOES-Next attitude estimation could be carried out using either 
a batch least-squares estimator or a Kalman filter. The Kalman filter is preferred 
because it is more straightforward to account for process noise with this method. 
The observability analysis of Section 4 shows that yaw data are needed to estimate 
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the constant yaw torque. Thus, a batch estimator must use at least one-third of 
an orbit, or 8 hours, of data to estimate the GOES-Next yaw across the period of 
yaw data outage. It is very likely that the dynamic models are not accurate enough 
to propagate across this interval without accounting for process noise. Thus, only 
the Kalman filter will be considered in this paper. 

The Kalman filter propagates estimates of the state x^ ,(+) and covariance matrix 
P k ^+) immediately after the (k-l) iX measurement to the time t^ of the k tn meas- 
urement by means of the following equations: 


X k <-> 


wt k- ‘k-l> %-!<♦> 



♦(t,., t') u(t') dt' 


(31) 


and 


p k ( -> W p k-i (+) * <V Vi) 

t,. 


(32) 


/ 


+ / $(t k , t') Q 1>‘(t k , t') dt' 


where the transition matrix $ is given by Equation (15) and the process noise 
spectral density matrix Q is defined by 

E[w(t) w T (t')] - 0 6(t - t') 


(33) 
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E denotes the expectation value, w is defined by Equation (12b), and S(t - t') 
denotes the Dirac delta (unit impulse) function. Equation (31) is simply Equa- 
tion (13) without the unknown process noise term, w. The explicit form for the 
process noise spectral density matrix is 


with 


and 




0 

V " 2 

o 

0 

0 


0 0 0 

0 0 0 

0 0 0 

0 <q z * V /h2 " q w /H 
0 -q w /H q w 


^4x4 " d ^ ag ^px* q pz’ q cx’ q cz^ 


(34a) 


(34b) 


(34c) 


where the latter notation means that Q 4x4 is a diagonal matrix with the indicated 
arguments as the elements on the main diagonal. The scalar spectral density q x is 
defined by 


ECw x (t) w x (t')] - q x 6(t - t') (35) 

and similar relations hold for q z , q w , q px , q pz , q cx , and q cz . The results of 
Section 3 allow closed form evaluation of the Integral in Equation (32), giving a 
very efficient means of covariance propagation. 
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When a measurement is processed, the state estimate and covariance matrix are 
updated as follows: 


M+) 


£.(-) + < k Cg k - G|/x u (-)] 


k k' 


(36) 


and 


P k (+) - (I g - K k G k ) P k (-) (I g - K k G k ) T + K k R k K^ (37) 

where g k is the vector of measured values at time t k , G k is the matrix relating the 
measurement to the state as in Section 4, R k is the measurement covariance and K k 
is the Kalman gain: 


K k- p k ( -> G J[w- )G J + R kj" < 38) 

If the attitude estimates are not needed in near-real time, an optimal filter- 
smoother may be preferable to a Kalman filter (Reference 12), but that option will 
not be considered in this paper. 

6. ACCURACY ESTIMATES 

In the absence of detailed simulations, quantitative estimates of the accuracy 
attainable with the proposed estimation procedure require the use of approximate 
models. This paper will consider the accuracy during periods when no yaw measure- 
ments are available, because these periods test the attitude estimation process 
more severely than periods containing yaw data. The observability analysis of 
Section 4 shows that the external torques are not completely observable without 
yaw data; the model will therefore be simplified here by eliminating the torques 
from the state vector. The remaining five-component state has an estimation error 
represented by a 5-by-5 covariance matrix, and the accuracy estimates will be ob- 
tained by computing an approximation to this matrix. The computation can be fur- 
ther simplified by averaging the covariance propagation Equation (32) over a 
nutation period, using the results of Section 3 for the transition matrix, and 
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then considering the limit at which the orbit rate Is negligible compared to the 
nutation rate. This permits the deletion of the roll rate and yaw rate from the 
state, resulting In an effective 3-by-3 covariance matrix. The steady-state co- 
variance can then be found by solving a quartlc equation, which Is possible In 
principle but inconvenient In practice. However, the simplified approach pre- 
sented below gives equivalent results for the steady-state covariance. 

The validity of the approximate steady-state covariance analysis rests on two 
quantitative aspects of the GOES-Next attitude determination. First, the orbit 

_3 

rate is much less than the nutation rate, by a factor of 2 x 10 . Second, the 

wheel tachometer measurements are so accurate that the yaw component of the inter- 
nal angular momentum, which is a component of the state vector, is essentially 
determined by tachometer measurements alone; the correlation between its errors 
and the attitude errors is thus effectively zero. 

Consider the dynamic Equations (7a) and (7b) in the limit that « n becomes infinite. 
The resulting equations will describe motion on time scales large compared to the 
nutation period, with nutation averaged out. The second derivative terms in Equa- 
tions (7a) and (7b) are seen to be negligible if nutation is ignored. The factor 
1/1 is replaced by « n /H using Equation (8), and H is held at its physical value as 
« n is taken to infinity. This means that the rotational inertia of the spacecraft 
body is dominated by the bias angular momentum on time scales larger than the nuta- 
tion period. The limits of Equations (7a) and (7b) are 

0 ■ « 0 r + y + (N)( + w 0 h)/H (39a) 

0 - « 0 y - r + (N z - N w )/H (39b) 

The two-component state 

x = [r + h/H, y] T (40) 

obeys the dynamic equation 

x-Fx+u+w (41a) 
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with 


and 



w - H 1 [w z , - w x ] T 


(41b) 


(42a) 


(42b) 


Equation (4) was used in deriving these equations. With the assumption that q x 
and q 2 , defined in Equation (35), are equal, the process noise spectral density 
matrix, defined by an equation analogous to Equation (33), is 

0 - (q/H 2 ) I 2 (43) 

where q denotes the common value of q x and q 2 - 

In this derivation, q is the spectral density of the process noise representing 

unmodeled external torques. The more rigorous method of averaging over a nutation 

period, as discussed at the beginning of this section, shows that a contribution 

1/2 q w from the momentum wheels should be added to q, the factor of 1/2 arising 
w 2 2 

from the time average of sin w n t and cos u n t. Since q w arises from torque rip- 
ple, as discussed in Section 2, its numerical value can be estimated as 


q w - 2(sin 1.656°) 2 (6 x 10 _7 /2ir)(N*m) 2 s - 1.3 x 10 _10 (N*m) 2 s (44) 

where the factor of 2 appears because two wheels contribute to h, and 2ir is the 
conversion from hertz to radians per second. The environmental torque errors are 
not well approximated by a white noise process, but the spectral density of these 
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errors can be estimated by multiplying the mean square torque errors by a correla- 
tion time. The maximum expected torque errors are about 10 percent of the ampli- 
tude of equinox torques shown in Figure 2. The quarter-orbit period, which also 
gives the filter memory span required to estimate yaw, is a reasonable estimate of 
the correlation time. Thus, an upper limit of the spectral density of external 
torques is 


q = ^ 5 x 10" 6 N«mJ (21541 s) = 5.4 x 10 -7 (N»m) 2 s (45) 

The wheel torque ripple spectral density, q w , is negligible compared to this. 

The speed of each momentum wheel is measured at 0.512-second intervals with error 
variance: 

R tach = [<0.0163) 2 /12 + (0.0978/3) 2 ] (2ir/60) 2 (rad/s) 2 = 1.2 x 10" 5 (rad/s) 2 (46) 

where the factor of 12 converts quantization error to variance and 2ir/60 converts 
revolutions per minute to radians per second. The roll is also measured at 
0.512-second intervals with error variance: 

R roll " + (0.01 ) 2 ](ir/180) 2 rad 2 = 3.3 x 10 -8 rad 2 (47) 

where ir/180 converts degrees to radians. These measurements can be combined to 
give a measurement of the first component of >T with error variance 

R k - "roll " 2 ('w sln 1-656 °/h) 2 R tac „ * R ro „ (48) 

since the tachometer errors are negligible. 
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Since the measurement interval At is much less than the nutation period, which 
has already been neglected in deriving Equation (41a), the measurements can be 
treated as continuous rather than discrete. This leads to a first-order differ- 
ential equation for the covariance (Reference 12): 


o* ryv ~ /%/vT _ 1 

F ■ + PF + Q - PG R GP 


(49) 


where 


G - [1 0] 


(50) 


and 


R = R k At = 1.7 x 10" 8 rad 2 s 


(51) 


Writing 


P = 


11 

P 12 

12 

P 22 


(52) 


and using Equations (41b) and (43) gives the three scalar equations: 


'll 


22 


. 2 « 0 p 12 + q/H 2 - p^/R 

(53a) 

-2 « 0 P-j 2 + q/H2 “ p 12 /R 

(53b) 

> 0 (p 22 “ P 11 5 ” p ll P 12 /R 

(53c) 
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Solving Equations (53a), (53b), and (53c) for the covariance for which the time 
derivatives vanish gives the following steady-state covariance: 


P 11 “ “ 0 RC(a + 3)(a - 1)]1/2 (54a) 

- 0) R a[(a + 3)(a - 1)] 1/2 (54b) 

r cz o 

Pl 2 = « Q R(a - 1) (54c) 


where 


a 


1 + q/ 



(54d) 


Taking a equal to the negative of the square root in Equation (54d) also gives a 
steady-state solution, but this is unacceptable because it gives a negative value 
for P 22 » which must be nonnegative. 

The principal quantity of interest for this analysis is the variance of the yaw 

2 2 5 

estimate, or p 22 - For GOES-Next, q/(o Q H R) = 3.9 x 10 » 1 , so to an excellent 
approximation: 

p 22 " “ 0 Ra2 ~ q/ (“o H2 ) (55) 


It is remarkable that the yaw variance is independent of the measurement error R 
in this limit. Equation (55) provides a general parameterization of the yaw ac- 
curacy as a function of the torque modeling errors. Inserting the GOES-Next 
values for the parameters gives 

c^w - <P 2 2 )1/2 ~ 6 ‘ 9 x 10 " 4 rad ■ °- 040 de 9 (56) 

This is very close to the performance expected by the spacecraft contractor (Ref- 
erence 10). 
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7. CONCLUSIONS 


Dynamic propagation promises to be a valuable complement to gyro propagation for 
GOES-Next, permitting observation of yaw during Sun sensor data gaps and estimation 
of environmental torques for setting the trim tab. Provided the attitude remains 
close to Earth pointing, the propagation equations are nearly linear. The added 
assumption that the body angular momentum is much smaller than that of the momentum 
wheels allows closed-form solution of those equations. The solution has terms that 
vary at the orbital rate and at the much higher nutation rate. Neglecting errors 
with greater than the nutation frequency allows closed-form expressions for the yaw 
accuracy, also. 

Error enters the propagation through imperfectly modeled environmental torques, 
control torques that are ignored because of lack of information, and random fluc- 
tuations in the wheel -bearing torque. The largest sources of error are expected 
to be the solar torque modeling and the neglected magnetic control torques. The 
stabilizing pitch momentum bias, the accurate wheel tachometer data, and the rel- 
atively small torques at geosynchronous altitude permit propagation in the pres- 
ence of these errors. 

Since the error propagation equations are linear, the Kalman filter proposed here 
is also linear and should be easier to develop than a general extended Kalman fil- 
ter. Hi th the transition matrices calculated in closed form, the filter should 
also be efficient to operate. In practice, this filter would process full orbits 
of data to estimate torque parameters and propagate the yaw during Sunless periods. 
The process and sensor noise levels expected for GOES-Next yield an estimated yaw 
accuracy of 0.040 deg following the gap in Sun coverage. This number is consistent 
with manufacturer estimates of yaw drift over that period. Successfully applying 
dynamic propagation to GOES-Next operations support would break new ground for the 
Flight Dynamics Division of the Goddard Space Flight Center and would provide a 
much needed backup to the usual gyro propagation. 
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